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A  three-dimensional  electrothermal  model  has  been  developed  to  investigate  the  spatial  resolution  of 
the  scanning  thermoelectric  microscopy  (SThEM).  We  found  that  if  the  electrical  resistivity  of  the 
sample  changes  abruptly,  the  SThEM  will  measure  a  voltage  close  to  the  local  thermoelectric 
voltage  where  electrical  resistivity  is  relatively  low,  rather  than  a  simple  weighted  average  of  the 
thermoelectric  voltage  distribution  based  on  the  temperature  profile.  This  is  due  to  the  presence  of 
internal  currents  in  the  sample.  The  spatial  resolution  of  the  Seebeck  profiling  is  limited  by  the  finite 
value  of  the  phonon  mean  free  path  of  the  sample  and  the  tip  size  of  the  microscopy.  With  a  tip  size 
around  1  nm,  the  scanning  thermoelectric  microscopy  can  achieve  a  spatial  resolution  of  the 
physical  limit  defined  by  the  statistical  nature  of  the  charge  carrier  and  phonon  behavior  in  a  very 
small  region.  ©  2005  American  Institute  of  Physics.  [DOI:  10.1063/1.2008381] 


Thermoelectric  devices  can  remove  heat  from  hot  elec¬ 
tronic  and  optoelectronic  chips  or  convert  it  to  electricity 
without  a  moving  part.  The  performance  of  a  thermoelectric 
device  is  mostly  determined  by  the  dimensionless  hgure  of 
merit  ZT.  ZT  is  defined  as  ZT={S^cr/ k)T,  where  S  is  the 
Seebeck  coefficient,  and  cr  and  k  are  the  electrical  and  ther¬ 
mal  conductivity,  respectively.  Recently  the  search  for  high 
ZT  materials  has  focused  on  nanoscale  structures  due  to  the 
enhanced  phonon  scattering  and  improved  thermoelectric 
power  factor.  ’  High  ZT  values  of  up  to  2.4  have  been 
achieved  for  the  thin  films  and  quantum  dot  superlattices. 

Since  the  physical  parameters  S,  cr,  and  k  are  modified 
locally  due  to  the  nanoscale  features  of  the  materials,  the 
conventional  characterization  methods  for  bulk  materials 
may  not  give  accurate  results.  Scanning  thermoelectric  mi¬ 
croscopy  (SThEM)  provides  a  way  to  measure  the  Seebeck 
coefficient  directly  with  nanometer  spatial  resolution.^  The 
setup  for  Seebeck  coefficient  prohling  is  shown  in  Fig.  1.  In 
an  ultrahigh-vacuum  environment,  the  scanning  tunneling 
microscope  (STM)  worked  in  the  constant  current  mode.  At 
each  point  of  the  measurement,  the  STM  current  loop  was 
disconnected  and  the  thermoelectric  voltage  measurement 
loop  was  connected.  The  STM  tip  stepped  in  by  the  piezo¬ 
controller  and  achieved  a  nanocontact  with  the  sample.  The 
sample  heater  caused  a  temperature  difference  of  30  K  be¬ 
tween  the  sample  substrate  and  the  STM  tip.  A  temperature 
gradient  was  developed  in  the  sample  and  localized  primarily 
in  the  nanocontact  region.  The  thermoelectric  voltage  V  was 
measured  and  the  local  Seebeck  coefficient  could  be  calcu¬ 
lated  by  S=V/AT. 

The  measured  Seebeck  profile  of  a  GaAs  p-n  junction 
showed  a  clear  Seebeck  coefficient  sign  change  at  the  tran¬ 
sition  from  the  p  region  to  the  n  region.^  However,  the  mag- 
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nitude  of  the  thermoelectric  voltage  at  the  transition  region 
was  much  lower  than  the  one  expected  from  one  dimensional 
model  based  on  the  measured  dopant  prohle.  This  is  attrib¬ 
uted  to  the  finite  sizes  of  the  tip-sample  nanocontact  and  the 
temperature  gradient  around  the  tip.  Because  heat  diffuses  in 
three  dimentions  from  the  sample  to  the  tip,  the  temperature 
gradient  and  Seebeck  voltage  are  distributed  in  the  sample, 
making  up  a  distributed  electrical  voltage  and  resistor  net¬ 
work.  The  measured  thermoelectric  voltage  is  the  terminal 
open  circuit  voltage  of  this  network.  In  this  paper,  we  present 
a  three-dimensional  (3D)  electrothermal  modeling  of  the 
Seebeck  profiling  of  a  GaAs  p-n  junction  by  scanning  ther¬ 
moelectric  microscopy  and  investigate  its  spatial  resolution 
for  the  first  time. 

The  thermoelectric  voltage  source  distribution  is  deter¬ 
mined  by  the  temperature  profile  and  local  Seebeck  coeffi¬ 
cient:  Eix,y ,z)  =  Six)AT{x,y ,z).  In  the  bulk  materials,  the 
heat  conduction  is  described  by  continuum  equation.  How¬ 
ever,  the  phonon  scattering  mechanism  must  be  taken  into 


FIG.  1.  Schematic  of  the  SThEM  setup. 
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account  for  nanostructures  because  the  phonon  mean  free 
path  is  at  about  the  same  scale  as  the  structure’s  dimensions. 
Thus,  we  treat  the  solution  of  the  temperature  profile  in  two 
regions  separately  in  this  model:  Around  the  tip-sample 
nanocontact  and  at  a  distance  of  a  phonon  mean  free  path  to 
the  contact,  we  assume  the  temperature  does  not  change 
since  there  is  very  little  phonon  scattering  in  such  a  short 
length.  Out  of  this  region,  the  temperature  change  is  gov¬ 
erned  by  continuum  equation  of  heat  conduction.  The  tem¬ 
perature  profile  also  depends  on  the  shape  of  the  nanocon¬ 
tact.  In  our  model,  we  assume  the  contact  is  a  round  and  flat 
plate  on  the  surface  of  the  sample.  In  a  phonon  mean-free- 
path  distance  to  the  plate,  the  temperature  is  assumed  to  be 
the  same  as  the  plate.  Out  of  this  region,  the  temperature  is 
solved  by  the  finite  difference  method  according  to  con¬ 
tinuum  heat  conduction.  The  as-obtained  temperature  distri¬ 
bution  is  a  first-order  approximation  of  the  actual  one,  which 
can  be  obtained  using  a  molecular  dynamics  or  Monte  Carlo 
simulation  method. 

The  Seebeck  coefficient  can  be  calculated  from  the  dop¬ 
ant  profile.  We  use  the  doping  density  from  the  experimental 
data  (Afl  =  9X  lO'*  cm"^,  N^=l.l  X  lO'®  cm“^).  The  electron 
and  hole  densities  can  be  solved  numerically  from  the  Pois¬ 
son  equation 


dx^ 


- [n(x)  -p{x)  +Na{x) 
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The  Seebeck  coefficients  from  the  contribution  of  electrons 
and  holes  are 
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where  F,{  rj)  is  the  Fermi-Dirac  integral  of  order  r.  For  tem¬ 
peratures  less  than  the  Debye  temperature  (344  K  for  GaAs), 
the  electron  and  hole  scattering  parameters  r^=ri,=  l/2.  The 
total  Seebeck  coefficient  including  the  contribution  of  both 
electrons  and  holes  is 


•S'(x): 


n{x)fi^Se{x)  +  p{x)fii,Shix) 

nix)p-e  +  pU)p-h 


(6) 


The  Schottky  junction  of  the  tip-sample  contact  changes 
the  band  bending  and  Seebeck  coefficient  profile.  However, 
for  a  nanocontact,  the  depletion  region  of  the  metal- 
semiconductor  junction  is  also  small.®  If  it  is  smaller  than  the 

phonon  mean  free  path,  it  will  only  have  trivial  effects  on  the 
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FIG.  2.  The  comparison  of  the  thermoelectric  voltage  calculated  based  on 
3D  simulations  for  different  phonon  mean  free  paths  and  the  ID  Seebeck 
profile. 


thermoelectric  voltage  measurement.  For  simplicity,  we  did 
not  include  the  effects  of  tip-sample  Schottky  junction. 

We  divided  the  planar  p-n  junction  into  48  X  48  X  24  cu¬ 
bic  cells  of  resistor  and  voltage  source  network.  The  resistor 
network  is  obtained  from  the  carrier  density  and  the  local 
conductivity  of  the  depletion  region  and  the  voltage  sources 
are  derived  from  the  3D  temperature  distribution  and  the 
local  Seebeck  coefficients.  A  system  of  linear  equations 
based  on  Kirchhoff’s  laws  (conservation  of  current  and  volt¬ 
age)  are  then  written,  their  solution  gives  the  node  potentials. 
From  this,  we  can  calculate  the  measured  thermoelectric 
voltage  at  the  STM  tip. 

For  a  highly  doped  (lO'®  cm“^)  GaAs  sample,  the  pho¬ 
non  mean  free  path  can  be  smaller  than  4  nm.  We  first  ne¬ 
glected  the  size  of  the  tip  contact  area  and  changed  the  pho¬ 
non  mean  free  path.  The  comparison  of  3D  thermoelectric 
voltage  simulation  for  different  phonon  mean  free  paths  and 
the  one-dimensional  (ID)  Seebeck  profile  is  shown  in  Fig.  2. 
The  magnitude  of  the  3D  curves  is  lower  than  the  ID  curve 
and  this  explains  why  the  measured  Seebeck  voltage  is 
lower  than  the  ID  simulation.  It  can  be  seen  that  the  magni¬ 
tude  of  the  3D  thermoelectric  voltage  profile  increases  with 
the  decrease  of  the  phonon  mean  free  path  and  converges  to 
the  ID  result.  This  indicates  materials  with  a  smaller  phonon 
mean  free  path  have  a  higher  spatial  resolution  of  Seebeck 
profiling. 

In  Fig.  3,  we  remove  the  effects  of  finite  phonon  mean 
free  path  by  setting  it  to  zero  and  change  the  tip  radius  alone 
from  1  to  4  nm.  It  can  be  seen  that,  as  expected,  a  smaller  tip 
size  gives  a  higher  spatial  resolution  of  Seebeck  profiling.  In 
Fig.  4,  we  combined  the  effects  of  finite  phonon  mean  free 
path  and  tip  size  and  set  the  phonon  mean  free  path  to  4  nm 
for  tip  radii  of  1  and  2  nm.  It  shows  the  spatial  resolution  is 
further  reduced.  The  measured  data  are  shown  as  dots,  which 
have  similar  magnitudes,  but  the  peak  location  is  shifted.  The 
reason  for  this  shift  is  not  very  clear.  It  could  be  due  to 
nonuniform  doping  distribution. 

It  is  interesting  to  note  that  the  measurement  of  the  local 
thermoelectric  voltage  is  not  a  simple  average  of  the  thermo¬ 
electric  voltage  in  the  region  covered  by  the  tip  size  and  the 

phonon  mean  free  path.  This  can  be  seen  from  the  fact  that 
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FIG.  3.  The  comparison  of  the  thermoelectric  voltage  calculated  based  on 
3D  simulations  for  different  tip  sizes  and  the  ID  Seebeck  profile. 

the  3D  simulation  results  are  lower  than  the  ID  Seebeck 
profile  almost  everywhere  in  the  device.  (The  3D  curves  and 
the  ID  curve  should  give  the  same  result  at  a  distance  far 
away  from  the  junction).  To  explain  this,  we  plot  the  varia¬ 
tion  of  the  local  thermoelectric  voltage  and  the  local  resis¬ 
tance  of  a  single  grid  cell  as  a  function  of  the  distance  from 
p-n  junction  in  Fig.  5.  When  we  have  a  temperature  distri¬ 
bution  that  extends  over  several  nanometers,  the  contribution 
of  local  thermoelectric  voltage  will  be  modified  due  to  the 
nonuniform  resistor  network.  As  a  simple  example  in  Fig.  5, 
if  we  have  a  voltage  source  Ei  in  series  with  a  resistor  Ri  at 
position  1  and  a  voltage  source  £2  in  series  with  a  resistor  R2 
at  position  2,  assuming  the  tip  contact  is  in  the  middle  of 
position  1  and  2,  the  parallel  combination  of  these  two 
sources  will  produce  the  following  open  circuit  voltage: 


FIG.  4.  The  3D  thermoelectric  voltage  simulation  for  the  finite  tip  size  and 
the  phonon  mean  free  path.  The  dots  are  experimental  data. 


Position  (nm) 


FIG.  5.  The  ID  thermoelectric  voltage  and  resistance  of  a  single  grid  cell  as 
a  function  of  distance  near  the  p-n  junction. 


Ri+  Ri 

This  is  due  to  the  internal  current  flow  in  the  structure. 
Since  the  magnitude  of  the  single  grid  resistance  decreases 
much  faster  than  that  of  the  ID  thermoelectric  voltage  as  we 
move  away  from  the  junction,  one  should  note  that  the  ter¬ 
minal  voltage  is  mainly  determined  by  the  lower  electrical 
voltage  £2  which  is  farther  from  the  junction  and  in  series 
with  the  smaller  resistance. 

In  summary,  using  a  3D  model  of  the  thermoelectric 
voltage  measurement  by  SThEM,  we  reveal  that  the  3D  ef¬ 
fects  cannot  be  simplified  as  an  average  of  the  ID  model  of 
thermoelectric  voltage  in  the  adjacent  region  of  the  tip- 
sample  contact.  The  measured  value  is  the  terminal  voltage 
of  a  network  of  thermoelectric  voltages  and  electrical  resist¬ 
ers.  The  spatial  resolution  of  this  nanoscale  measurement 
technique  depends  on  the  phonon  mean  free  path,  and  the 
size  of  the  tip-sample  contact.  One  should  note  that  the  local 
thermoelectric  voltage  can  only  be  defined  and  calculated 
within  a  uniform  volume  larger  than  the  relaxation  lengths  of 
the  energy/charge  carriers  and  the  inter-dopant  spacing  due 
to  its  statistical  nature.  Thus,  this  nanoscale  Seebeck  mea¬ 
surement  technique  can  almost  achieve  the  spatial  resolution 
of  the  physical  limit  (the  contact  size  can  be  made  as  small  as 
1  nm  and  below). 
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